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While the possible roles of GRBs' progenitors as Gravitational Waves (GW) sources have been 
extensively studied, little attention has been given to the GRB jet itself as a GW source. We expect 
the acceleration of the jet to produce a Gravitational Wave Memory signal. While all relativistic 
jet models display anti-beaming of GW radiation away from the jet axis, thus radiating away 
from directions of GRBs' 7 radiation, this effect is not overwhelming. The decrease of the signal 
amplitude towards the cone of 7-ray detectability is weak, and for some models and parameters the 
GW signal reaches a significant amplitude for much of the 7-ray cone. Thus both signals may be 
jointly detected. We find different waveforms and fourier signatures for uniform jets and structured 
jet models - thus offering a method of using GW signatures to probe the internal structure and 
acceleration of GRB jets. The GW signal peaks just outside the jet (core) of a uniform (structuted) 
jet. Within the jet (core) the GW signal displays wiggles, due to a polarization effect; thus for 
a uniform jet, the peak amplitude accompanies a smoother signal than the peak of a structured 
jet. For the most probable detection angle and for typical GRB parameters, we expect frequencies 
< GOOHz and amplitudes h ^ 10~^^. Our estimates of the expected signals suggest that the signals 
are not strong enough for a single cluster of DECIGO nor for aLIGO's sensitivities. However, 
DECIGO's 2-clusters mode should detect such signals from most detected GRBs, as well as from 
events whose GRB's are beamed away from us (including orphan afterglows). 



I. INTRODUCTION 

General Relativity predicts Gravitaional Waves ([T), 
and evidence of their effect on astrophysical systems has 
established their existance ([IHl])- However, no direct 
observation of such waves has yet been achieved. The 
searches for gravitational wave sources ([5]) can be sim- 
plified if one knows in advance and with enough precision 
where and when to look. When searching for a GW sig- 
nal amidst background noise, it is also helpful to know 
external parameters regarding the source that might help 
characterize it and its wave form- template. In addition, 
the ability to study a single astrophysical event simulta- 
neously both in the EM spectrum and in its GW spec- 
trum can provide much more information and insights 
than the mere sum of the separate measurements. For 
these reasons, it is interesting to look for Astrophysical 
objects that emit both EM and GW radiation, and it 
is interesting to explore whether and under what condi- 
tions they are jointly detectable (|6H17j). GRBs are the 
most natural candidates observed to emit strong EM ra- 
diation and expected to emit GW radiation as well [7]. 
In particular, short GRBs are expected to be associated 
with merging binary neutron stars, which are the classi- 
cal candidates for detection by advanced GW detectors 

Different proposed mechanisms concerning GRB's have 
been suggested and studied as possible sources for Grav- 
itational Waves. However, little attention has been given 
to the relativistic jet itself as a source. Regardless of 
the progenitor type, all GRB's share an initial phase of 
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an accelerating relativistic jet. GRB's involve explosions 
which release about Ej ^ lO^^ergs of energy in acceler- 
ating jets (a forward jet and a backward jet), reaching 
ultra-relativistic velocities with Lorentz factors F >100 , 
or even up to 1000 [23 . The jets are strongly beamed 
into narrow (double) cones, of angle ^0 0.1. The dura- 
tion of the observed bursts ranges from less than O.Olsec 
to more than lOOsec, while their light curves show rapid 
variability, at times on scales of about O.lsec. The ultra- 
relativistic non spherically-symmetrical acceleration of 
this energetic jet is expected to produce Gravitational 
Radiation ( [24H33] ). In this work, we focus on this GW 
signal, from a single pulse of an accelerating jet, regard- 
less of the progenitor. 

We model the jet in 3 stages of increasing complexity: 
first, as an instantaneous infinitely narrow jet pulse (the 
"point particle approximation", section (2)) that enables 
us to estimate the signal strength; second, as an instan- 
taneous cone with a finite width and mass distribution 
(aggregate models, section (3)) for an understanding of 
the angular distribution; and finally as a cone of both 
finite width and finite time (a prolonged acceletation 
model, section (4)) for am estimation of the temporal 
structure (waveforms). We analyze the results of each 
model, determining the directions of radiation (and in 
particular the possibility of a joint detection of GW and 
EM (7-ray)), and characteristic signal amplitudes, en- 
ergy output, and frequency behaviour. We build each 
model using results from the simpler ones. In section (5) 
we discuss the possibility of detection, and future work. 
Note on conventions: Adhering the General Rela- 
tivists' convention of "natural units", we use (c=G=l). 
However, certain formulae for comparison with physical 
results and scales employ the cgs system. The space time 
metric signature is diag{ — h ++)• 
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II. INSTANT ACCELERATION OF A 
POINT-PARTICLE 

A. GW Memory and the Point Particles 
Approximation 



This is a purely quadrupole radiation pattern, analo- 
gously to the Lienard-Wiechert Potentials for electro- 
magnetic radiation [39] . 



While the progenitors emit quasi-periodical gravita- 
tional waves (Coalescing of Compact Binaries, etc.), the 
acceleration of the jet itself is prompt and non-periodic, 
and thus we do not expect "waves". Rather, we expect 
the jet to produce a gravitational wave memory (GWM; 
or Zero Frequency Limit, ZFL [25H34] ). defined as: 



{t)-\im /i+,x(t), (III) 



where /i+ and hx denote the plus and cross standard 
polarizations of the wave function [35) . 
Einstein's equation reduces - in the regime of linearized 
gravity and under the harmonic gauge condition [36 - 
to a wave eqaution for the metric perturbation (/i^^^ = 
9fiu — Vfiu)^ with the Stress-Energy tensor serving as the 
source term: 



1 



-167rV = -167r(T^, - (11.2) 



The solution of this equation is: 



(11.3) 



The simplest approximation of a Stress-Energy tensor 
for a GWM calculation is obtained using the "Collision 
Approximation" , which models the process as an instan- 
taneous exchange of ingoing point-particles with outgo- 
ing point-particles. The gravitational signal is derived 
from the incoming/outgoing momenta of the n particles 
(P^/P^^) giving the stress-energy tensor T^^: 



n E„ 



(11.4) 



Combining these two equations and integrating over 
the ^-functions gives the gravitational signal. Simi- 
lar "Collision Approximation" methods have been used 
by Weinberg [36 , Piran[9 and others (j26]-[29],07],g8j). 
Smarr [34 has shown that this form reproduces (using 
momentum conservation) the gravitational wave signal 
form used by Braginsky, Thorne and others for instanta- 
neous collisions (of N particles with rest-masses rriA and 
velocities va at angles 9a and distance r to the observer, 
and with the A marking the memory effect as defined in 
eq. dint): 
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B. The Jet as a Point-Particle 

The simplest approximation of the jet's GW signal 
would be to treat the entire jet as a single point-mass 
undergoing an instantaneous linear acceleration from rest 
to r, reaching the total jet energy Ej. This model is of 
course over-simplified, as it neglects both the jet's finite 
width (the model assumes it to be infinitely concentrated 
and narrow) and the finite duration of the acceleration 
phase (which the model assumes to be instantaneous). 
However, it provides an understanding of the relevant 
characteristic signal and output energy scales (computed 
from the accelerated jet energy Ej)^ as well as an approx- 
imate angular distribution. 

The change induced in the metric between t ±oo for 
an observer at a distance r away is determined by the an- 
gles between his line-of-sight and the various velocities, 
the distance and the energy of the jet. As in [26l [27] we 
set N=2 for an initial "incoming particle" at rest [Pi = 0) 
and a final "outgoing particle" of mass m/ = m = Ej /T 
with = /3(sin cos ^, sin sin (/), cos ^), where P and the 
angles are in the observer frame (the z axis points to the 
observer). The contribution of a point-particle acceler- 
ating instantaneously from to /3 to the gravitational 
perturbation h is 



AhPP = Ah^ ^ iAh, 



2Vmp^ s'm^e ^^^^ 



1- pcosO 



e''^. (n.6) 



We note here that for relativistic velocities {(3^1) the 
typical scale for the amplitude of the wave will be deter- 
mined by 



(n.7) 



(n.5) 



We see that the angular distribution of radiation from 
an instantaneously accelerating point source is set by 
two competing effects: the multipole pattern of the ra- 
diation, expressed as sin^ 6 (for multipole order n), di- 
rects the radiation away from the forward direction, 
while the transformation from the CMRF to the observer 
frame Lorentz-contracts the radiation forward, by a fac- 
tor of (1 — f3 cos 0)~^. For the GW quadrupole radiation 
{how ^ sin^ 0/{l — PcosO)) the result is "anti-beaming" 
away from the forward direction (fig. Q), i.e. an almost 
uniform distribution outside a forward beam of width 
e - (fQ^ ^^^^^ I _ i^^^^Q _ p-i) j24]^ This 

implies detection of both EM and GW radiation requires 
attributing a finite width to the jet. 
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FIG. 1. Antenna Diagram showing relative beamed direc- 
tional intensity (energy) for isotropic (red) and quadrupole 
(blue, dotted) radiation from a point source (infinitesimal 
width), for F = 7, /3 = 0.99. The quadrupole radiation pat- 
tern is the expected for GW radiation from a point source. 
The diagram also shows the relative directional signal ampli- 
tude of the GW signal (blue, dash-dotted). 



2. Physical Parameters 

Assuming a physical time-scale for the acceleration 
process, ST ^ SOmse^ we expect a waveform display- 
ing typical frequencies / ~ 30Hz in angles away from 
the axis of acceleration. As the angle is closer to the 
acceleration axis, the typical frequencies are larger, in- 
versely to the time contraction (eq. (ILIO)). This places 
these sources in the frequency band for detection with 
the earth-based detectors (such as LIGO, Virgo and for 
smaller frequencies - DECIGO and NGO). 
We consider characteristic jets reaching total kinetic en- 
ergy ~ 10^^ ergs and a final Lorenz factor F ~ 100. 
This sets the nor maliz at ion c ondition for the mass pa- 
rameter from (eq. (ITG), ( |IL7| )). Integrating over all solid 
angl es an d restoring the physical units and parameters 
(eq. (|II.7[)), the total output is: 



(11.12) 



1. The energy flux 

The energy carried by the GW depends on the pertur- 
bation amplitude as well as on the frequency. We have 
calculated the GW Memory (ZFL) effect of a GRB jet, 
considering only initial (rest) and final (accelerated) par- 
ticle momenta, or equivalently, assuming instantaneous 
acceleration. Thus for energy calculations we attribute 
a finite timescale for the acceleration 5T in the source 
frame. In the observer frame. 



(A/.)^ 



(II., 



and by integrating over observed frequencies up to -^^tjq^ , 
we find the radiated GW energy 



dn Stt ST' (9) ' 



(II.9) 



where the observed timescales are related both to the 
acceleration timescales and to the observation angle by 



5T'{0)=5T{l-/3cos0), 



(11.10) 



as pointed out by Segalis & Ori ^25j. Thus the angular 
dependence of the energy flux is (suppressing the prefac- 
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and factors of G and c) 
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Using 0(1) estimates for Ah^^^{Oy) outside the narrow 
forward cone, and substituting the physical parameters 
(for r ~ 500Mpc), we estimate the energy and signal 
amplitude as: 



T^PP 



i^pp 
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lO^^erg 
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lO^^er^ 
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M 
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(11.13) 



(11.14) 



III. INSTANTANEOUS ACCELERATION OF 
AN AGGREGATE SOURCE 

A. Aggregate Source Models 

We turn now to study a slightly more realistic model: 
an aggregate model of an expanding shell (a sector of an 
envelope, fig. ([2|). The ejected mass is now distributed in 
a shell according to some angular distribution f{6) that 
is axis-symmetrically around the jet's central axis. The 
jet axis is at an angle Oy relative to the observer's line- 
of-sight (the observer is at a distance r away). We treat 
every mass-element in the shell as an accelerated point- 
particle, and find its contribution Ah^^ to the signal re- 
ceived at r. For every Oy we integrate the contributions 
Ah^P over (^, (/)), to find the total signal from the jet, h^: 

Ah^{Oy) = J f{0,(l))AhPP{0,(l),Oy)smOdOd^, (IILl) 



The additional (1 — (3 cosO) term in the denominator en- 
hances the beaming of the energy flux in the forward di- 
rection, and its peak is still outside a 6> ~ a/2F~^ forward 
cone (fig. ([T])). 



^ we assume a stage of outward baryonic expansion from R ~ 
10^ cm to R ^ 10^ cm, where for most of the expansion the ve- 
locity V ^ c. 
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where / is a normalized for conservation of total ejected 
energy between the different models and parameters. We 
consider two possible models for the mass distribution in 
the jet (Both models are of course idealized approxima- 
tions to the real, unknown angular distribution of a GRB 
jet): 

a. A uniform jet: 



m 



{I 



e<eo, 
e>eo. 



b. A structured jet (following [40 ): 

1 9 < 9 rare = F ^, 
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FIG. 2. Our Jet modelled as a sector of an envelope of mass- 
energy, expanding. 



FIG. 3. Antenna diagrams showing angular distribution 
of the GW signal amplitude emitted from F = 100 for 
uniform and structured jets of different jet opening angles 
^0(0.01, 0.05, 0.1, 0.2). The axes mark relative directional am- 
plitude, the X-axis marks the jet axis and the y axis is perpen- 
dicular to it. The blue curve describes the uniform jet, while 
the red describes the structured jet. Each pair is normalized 
to the same total jet energy, demonstrating the relative wave 
intensities and angular distributions. Both models coincide 
for Oo = 0.01; for larger ^0, the structured jets are narrower 
and their radiation is both stronger and it is beamed into a 
sharper emission cone. The green line demarcates the opening 
angle of the jet itself, approximating the region of detect abil- 
ity of the 7-ray signal from the GRB jet. Small angles are 
zoomed-in (boxed), showing the regions of joint detectability. 



B. Results 



We have numerically calculated the integral eq. (III.l) 
for both jet models with the same total jet kinetic energy 
of Ej = 10^^ ergs. Fig. ^ depicts a graphical compar- 
ison betwenn the angular distribution of the GW sig- 
nal for uniform and structured jets, for typical param- 
eters (F = 100; 6>o = 0.01,0.05,0.1,0.2). Fig. ^) and 
^p) depict together the energy and angular distributions 
for uniform and structured jets of different ^0 (respec- 
tively). Note that the minimal openning angle consid- 
ered 0.01 corresponds to the narrowest possible jet, for 
which 00 = T-^. 

Combining our numerical results for Ah^^*{Oy) with 
all of th e physical parameters (for r ~ 500Mpc, using 
eq. ( |II.7| )), we find th e jet w idth's effects on the emitted 
signal (compare with ( II.14| )). We find that for a uniform 
jet: 



^peak "^2-10 



M 



lO^^erg J \0m 
while for a structured jet: 



-0.1 



500Mj9c 



(III.2) 
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500Mpc, 
(III.3) 



For the common parameter of F = 100, our models for 
a uniform and a structured jet coincide for Oq = Score = 
F~^ = 0.01. However, the signal strength and energy de- 
crease as Oq increases, with the uniform wave decreasing 
faster. This result is also evident in fig. ([3|,(|4|. This re- 
sult is general: for any Oq > 6 core = T"^, we expect the 
wave pattern from a structured jet to be stronger than 
the one from a uniform jet with the same paramers, and 
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FIG. 4. Antenna diagrams for the signal amplitude of GW 
emitted from uniform (left) and structured (right) jets of dif- 
ferent opening angles • All are normalized to the same total 
jet energy lO^^er^s. The axes mark relative directional am- 
plitude, the X-axis marks the jet axis and the y axis is perpen- 
dicular to it. The radiation pattern becomes wider, weaker, 
and more anti-beamned as increases. Both models coincide 
for 6*0 = 0.01. 
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FIG. 5. The region of a possible joint detection of a GW 
signal with 7-ray signal, relative to the total 7-ray detection 
cone. The left figure describes uniform jets, structured jets 
are on the right. For each F value the graph shows the results 
for 3 values of the threshold signal amplitude (rel. to the peak 
amplitude): 50% (bottom, dotted), 30% (middle, dashed) and 
10% (top, continuous). Different F are color-coded as red 
(F = 50), green (F = 100), blue (F = 500) and purple (F = 
1000) 



we expect the GW from uniform jets to decrease faster 
with ^0 than the corresponding structured ones. This is 
due to the higher spherical-symmetric character of the 
uniform jet: for 60 > F"^, the structured jet is always 
more focused than its uniform counterpart, and as ^0 in- 
creases, it remains more focused than the uniform. 
We defer more thorough study of possible detection by 
currently planned GW detectors (advanced LIGO, NGO, 
DECIGO, etc.), but do note that these signal amplitudes 
are weak (advanced LIGO [46] and DECIGO gllllH] of- 
fer strain sensitivities of ^ 10~^^ — 10~^^). Thus, the 
exact parameters of the jet, as well as the jet's internal 
structure model, might affect the possibility of detection. 



C. The Angular Distribution and Possible Joint 
Detection 



The GW signal from a single point-source is anti- 
beamed away from the axis of acceleration (at angles 
^peak ^ r~^/^). The relation between the detection cone 
for the GRB's 7-ray signal and the GW detection range 
for the uniform an d st ructured aggeragte models is more 
complicated (fig. (3j4ij4]3)). For a very narrow uniform 
j^t (^0 ^ r~^), the angular spread is similar to that of 



a point source, and Opeak ^ F~^/^. The peak angle in- 
creases (away from the jet axis) as 60 increases, and for 
wider unifrom jets we find approximately O^J^^ ex ^q"^^. 
Thus for all ^0 the uniform jet GW signal peaks above 



both Opeak 



> F~^/^ and 6>q"^^^. This is a consequence of 
anti-beaming of the GW radiation away from the source 
axis. As the 7-rays are only visible up to an observer 
angle of ^0 + F"^, which is smaller than Opeak ^ the peak 
of the GW radiation is generally outside of the 7-ray de- 
tection cone. The beaming effect is more pronounced at 
higher values of F; but it is qualitatively the same. Thus 
the GW peak cannot be observed by an observer that 
detects the 7-ray signal. 

However, the fall-off of the GW signal towards the jet axis 
(and the cone where its 7-rays are visible) is slow, as each 
of the point-particles' quadrupole signal itself is broad. 



Therefore if we are within the jet's 7-ray cone, we might 
still see the corresponding GW signal at less-than-peak 
amplitude. For a uniform jet detecting the peak 7-ray 
signal is possible anywhere within the 7-ray cone (up to 
maximum angle Oq + F~^ from the axis). The range of 
possible joint detectability, given as relative area (area 
of joint detectability region / total area of 7-ray cone) is 
shown in fig. ([5|. 

The situation is very different for structured jets. While 
we still consider their 7-ray signal to be as detectable 
anywhere up to ^0 + F"^, the GW signal comes mostly 
from the F~^ core, thus its peak is closer to the jet axis 
i^plak^^ oc 6>^-^^) and the signal fall-off away from the 
peak is slower. While in both models the radiation pat- 
tern becomes wider, weaker, and more anti-beamned as 
^0 increases, we find that this change is less pronounced 
for structured jets (fig. Q). As fig. ([5| shows, this im- 
plies a higher possibility of detection for structured jets, 
for the same parameter space. 

Having examined ranges of F values (10.. 1000), Oq values 
(0.01. .0.5), and relative threshold amplitudes (10%, 30%, 
50%) we find that even for uniform jets the result is quite 
optimistic, allowing joint detection for much of the 7-ray 
cone area (approximately 80-90%, 50-80%, 20-50% cor- 
responding to the 10%,30%,50% thresholds). Naturally 
the lower the threshold, the greater is the probability for 
joint detection. But even for a high 50% threshold (for 
which we loose only a factor of 2 in the detection rate) 
we still expect a considerable joint detection probabil- 
ity. We note also the dependence of the joint detection 
probability on F and on ^o- for smaller relativistic F fac- 
tors, beaming/anti-beaming effects are less pronounced, 
and the possibility of a joint detection is larger. Re- 
garding the opening angle Oq, we find the largest prob- 
abilities for a joint detection for very small values of ^0 
(approaching point sources; and also implying stronger 
GRBs). The probability decreases as ^0 increases, up to 
an angle ~ F~^/^, and then it rises very slowly rises. For 
structured jets, the parameter range of joint detection is 
even wider, allowing detection even of signals at the 50% 
threshold at over 80% of the 7-ray cone for reasonable 
models (^0 ^ 0-1) 7 and over a larger fraction of the range 
for larger Oq^ as shown in fig. ([5|. 

For both models, the availability of an EM trigger (7-ray) 
setting the time and position in the sky of the source can 
significantly increases the sensitivity of the GW detectpr 
as compared to a random search [7], thus increasing the 
chances of a joint detection. 
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IV. A PROLONGED ACCELERATION MODEL 

A. The Model 

So far we have considered jet models of different angu- 
lar distributions, undergoing instantaneous acceleration 
from rest (/3 = 0, T = 1) to /3 - c, T - 100 (following 
IS US [261 EH [38] ) . We used the duration of the accelera- 
tion ST just to estimate the typical frequencies. However, 
the instantaneous model is of course nonphysical, and in 
order to calculate waveforms, as well as the GW spec- 
trum, we must examine more specific models for the ac- 
celeration, where higher multipoles than the quadrupole 
appear. 

Our treatment of the GW radiation is analogous to the 
treatment of EM radiation in terms of Lienard-Wiechert 
Potentials (as in [39^), where the radiation contribution 
of an instantaneously accelerating point-particle source 
(with all quantities evaluated at retarded time) is given 

by 



V) X a] 



e 1 n X ((n 
c R {1 — V ■ n) 



(IV.l) 



In cases of linear acceleration such as ours, vxa = 0, and 
so n X (n X a) produces a sin^ term. Therefore for an 
instantaneous acceleration we retain only dipole EM ra- 
diation, and analogously only quadrupole GW radiation, 
relative to the instantaneous position at retarded time. 
When the acceleration is prolonged rather than instanta- 
neous, the position of the radiation source changes during 
the emission process. Therefore the radiation signal is 
a superposition of lowest-multipole terms from different 
times and sources - and thus includes higher multipoles. 
It is these contributions that we do consideJH] here. In 
general, the calculation of the radiation generated by a 
matter distribution p(t',x') and seen by an observer at 
(t, x) involves an integral of the form 



i?,x)=ff.»x'.(,-,'-|x-x'|).(,',xOA„-,x;x'). 

(IV.2) 

where Ah^^ is given by eq. (II. 6) for r = x — x', and the 



'dt 



delta function enables us to calculate using the retarded 
time t^e^(x, x') = t — |x — x'|: 

- {t, x) = /// d3xV(C x') A/i^^(x; x'). (IV.3) 



Inspired by the Fireball model [8l[4T], we expect the jet 
to accelerate from = 1 to F/ ^ 100, approximately 
linearl}j^ 



F (X r{t) - t. 



(IV.4) 



We use the total acceleration time Tf ^ 30ms. For sim- 
plicity we assume acceleration to be angularly uniform 
over the jet {0 < Oq): 



1 + 



< t < T/, a 



f 



■f 



T 



. - Tf 

(IV.5) 

At every instant t, we assume that all the ejected 
mass/energy is at the same radius r{t) from the center 
of the GRB, with the entire envelope- sector (fig. Q r] 
expanding radially at the same (angularly uniformj 
instantaneous velocity Pit) and Lorentz factor F(t). 
This model approximates the Radiative Fireball model 
[41 , where F cx r ~ t. This allows us to simplify 
eq. (IV. 2) by requiring 



(IV.6) 



using the jet-model functions / as in eq. (III.l). Thus 



dh ^'^^ ff 

— {t, ^)=JJdcos0'd^'f{0', 4>')rlet{(i', 4>', 0, t)AhPP{0', 4>', 0), 

(IV.7) 

where now the integration is performed over all the mat- 
ter elements responsible for the signal arriving at x at 
time with rret{0' ^ 0', 6>, t) describing the distance of each 
source element from the origin. 



Rat her th an calculating eq. (IV.7) by explicitly solving 
eq. (IV.6) and subtituting rret(0^(t)' ^O^t) in the integral, 
we have taken a different numerical route based on our 
earlier calculations. We divide the source in space over 
an angular grid (5000 x 100) of point-sources, and fol- 
lowed each such source as it accelerates (over 50000 time 
steps). Each point-mass accelerates linearly, maintaining 
its angular position (^,0), but at e ach tim e step its r,/3 
and F increase (F(t) linearly by eq. ( IV.5[ ), j3(t) matches 
F(t), and r{t) as a time integral over /3(t)). For each such 
point-mass at each time step, we model its acceleration 
as the disapperance of an incoming slow particle, and 
the appearance of an outgoing faster one. We sum both 
signal contributions using eq. ( |II.5| ), the positive (outgo- 
ing) contribution and the negative (incoming) one. We 



^ Higher multipoles can also appear if the acceleration is non- 
linear, i.e if the direction of a is different than the direction of 
V. Then the numerator aquires a new term— ri x {(y) x a), which 
is different from the lowest multipole term. Such terms are ir- 
relevant for our jet models, where the particles accelerate along 
straight lines. The method eq. ( |IV.1| is used to calculate the EM 
Larmor formula ( 39 ), including only the instantaneous dipole 
term, and the GW Quadrupole Formula(^2^). 



^ This is of course an approximation, and might break down in a 
realistic model, for example when the jet is propagating within a 
Collapsar or if the acceleration is magnetic rather than thermally 
driven. 

^ or rather partial envelope, for a uniform jet. For a non-uniform 
jet, such as the structured jet we examine, the mass density varies 
over the angle from the jet axis, but the envelope-front still all 
shares the same r{t) 
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record the net signal along with its time of arrival at the 
observer. The final observed signal at time t is the sum 
of the contributions arriving at that instant. 
Since Ah (xT oct (eq. (11.6) and eq. (IV.5)), each consec- 



utive contribution of the envelope to the signal increases 
with time (until the jet reaches F). The GW signal is 
"anti-beamed", outside an angle ~ r~^/^ away from 
the direction of acceleration, with F ex t, the signal emit- 
ted at t is anti-beamed outside an angle 6 ^ In 
the forward directions, where 7-ray radiation can be de- 
tected, we find that the GW signal contribution is larger 
at later times (when F~^/^ is small enough as not to 
encompass the observer), while the siganl contributions 
emitted earlier (when the observer's angle was ^ F~^/^) 
is undetectable. An observer in the forward (7-ray) cone 
would therefore miss the early signal, and only catch the 
later parts. The reverse-jet also contributes to the late 
signal (fig.([7]3)), as a shallow increase arriving later (un- 
til ~ 2T/'). However this signal is much weaker than the 
forward one (~ 10~^ compared with the forward signal, 
due to the beaming factor (1/(1 — /3cos6)). 
We expect to see the signal from the same Equal Arrival 
Time Surfaces (EATS) studied for EM radiation ([42]- 
|44J). These are characterized by two time-contraction ef- 
fects: radial time contraction (forward) and angular time 
contaction (sideways) [9 . The emitting jet travels at al- 
most the speed of light (in the observer's frame), so in 
the forward line of sight (FLOS) the jet and the signals 
it emits almost catch up with the signals emitted earlier. 
Thus the time delay between observing signals emerging 
at times ^1,^2 is contracted to 



At 



1 

2f2 



(^2-^1). 



(IV. 



Similarly (as for 7-ray waveforms [9^), the parts of the 
jet expanding at an angle off the FLOS are more distant 
from the observer relative to the parts of the jet along 
the FLOS. This introduces a time-delay that depends on 
the angle of the source (increasing away from the FLOS), 
and increases linearly with t, as ct is approximately the 
instantaneous radius of the envelope: 



At - {t2 - ti){l - cosO) - —{t2 - ti) (IV.9) 

in the arrival times between the signals from two sources 
emitting the signal simultaneously, on and of the FLOS. 
These two effects can be geometrically expressed together 



(a) 



Loieiitz Gainma factor r=100^ Uniform Jet Qo=0.1 

Observer 9^=0 01,0.05, 0.1,0 2,0.3, Q.5 




FIG. 6. h(t) for different observeration angles Ov (color- 
coded) ;r = 100, Oo = 0.1. Time is measured in dimensionless 
units t/T/, and h is scale d to h/ j^), where ^ - 10"^^ 
(compare eq. (II. 6), (II. 7), (II. 14)). (a) The entire waveform 
reaches plateau extending to 2T/ ^ 60ms. (b) A close-up 
of the high-rising signal from the forward jet, around r~^T/. 
(c) The very early signal, visible sooner away from the line of 
acceleration. 



Loreiitz Gamma factor r=100. Uniform Jet 8(5=0.1 

Observer G„=0.01. 0.05, 0.1,0.2,0.3, 0.5 




FIG. 7. /i as a function of time for different observeration 
angles ^^;;r = 100, = 0.1. (a) The waveform on a loga- 
rithmic time scale. The peaks are visible around 0.5r~^T/. 
(b) A close-up of the early signal, which is visible sooner the 
further off from the line of acceleration, (c) The late signal 
{t > l.STf) continues to rise very moderately, due to the back- 
wards jet. The final drop in h marks to the termination of 
the jet expansion, at 2T/. 



The Waveform 



+ F- 



-{t2-ti). 



(IV. 10) 



However, unlike the case of EM radiation, contributions 
to the GW signal from the FLOS or close to it are 
undetectable due to "anti-beaming". Thus > F~^/^ 
for any contribution detectable by the observer, and the 
contraction factor is always 0(F~^) rather than 0(F~^). 



We have run simulations of a uniform jet with a con- 
stant, uniform linear acceleration, over a wide parameter 
space of F, ^0, We have used different numbers of in- 
tegrations steps to confirm the numerical validity of the 
results. The resulting waveforms demonstrate the ex- 
pected effects. Tf factors out linearly, and the waveform 
results are identical (up to scaling) for bursts of differ- 
ent durations which reach the same final F; prolonging 
or shortening the acceleration is equivalent to stretch- 
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FIG. 8. GW signal from uniform jets of 6*0 = 0.1, T = 100, for 
different observation angles Ov •, around 0.1. (a) shows h{t) and 
(b) shows h{t). For Ov < ^o, shortly after the signal buildup 
begins h dips, reaching negative values. For sufficiently low 
6>^, h itself becomes negative (6*0 = 0.06,0.05,0.04,0.03). 
Later, the signal rises again. For Ov ^ Oq phase- variations 
causes wiggles to appear around the peak (seen for Oq = 
0.08,0.09,0.1,0.11) 



Comparison of Lorentz Gamma factor 
r=500,400, , . . ,150,100,75,50, 

Uniform Jet 80=0.1 

Obs Angle 9v=0.1 




FIG. 9. h(t) (a) and h ((b) - linear scale, (c) - log scale) for 
different Lorentz F values (500,400,300,200,150,100,75,50,25) 
(^0 = 0.1, Ov = 0.1). The signal is stronger, and peaks sooner, 
the higher the Lorentz factor, h is scaled to /i/(^^), and h 

to ii/im 



ing the waveforms in time - or shifting them in fourier 
space. We describe here t he results for a characteristic 
Tf = 30m5. Figures ( 6|7|8 ) show the waveform from dif- 
ferent observation angles for a uniform jet of parameters 
r = 100,6>o = 0.1. 

Fig.(|6|i) depicts the double- feature of the aggregated sig- 




typical time scales in h(t) 

Uniform Jet 6c.=0.1 
Lorentz Gaimua factor r= 100 



— first rise to peai</2 
— peal< 

last fa ii to peai</2 



FIG. 10. A comparison of timescales for ^ from a uniform 
jet of ^0 = 0.1, F = 100 as a function of the observation angle 
6v. Shown are the times of the peak of as well as the 
times when it reaches 1/2 of its peak value, both on the rising 
and on the falling edge. 



nal, from the forward and reverse jets. First, the con- 
tracted forward jet's signal (fig.([9fc)) rises rapidly over 



a time scale tw.sp ~ T ^Tj, matching our anticipation 



(eq. (IV. 10) and on). This is most clearly visible for small 
observation angles {Oy ^ 6>o), where the jet width sets the 
angular time-scale; for larger angles, the peak and total 
duration timescales are set by Oy > Oq (fig. ([lO|). It is 
followed by the long (~ 2Tf ~ 60ms) and weaker (by 
~ 10~^) backward-jet-signal (fig.([7]3)). Fig.([7]3) demon- 
strates clearly the anti-beaming effect and its dependance 
on time (via T{t)). The onset of the signal for the dif- 
ferent observers is determined by the observation angle 
Oyi The most off- axis observers observe the signal first. 
As r increases with relativistic beaming wins over, and 
observers closer to the forward axis detect the signal. 
Fig.([7]3) depicts the angular spreading of the forward 
signal. We find that for our characteristic Tf = 30ms 
burst, the strong signal from the forward jet lasts for 
~ 0.2 — 0.3ms. This angular spreading timescale (~ 
0.5(^o)^^/) is longer than the timescale of the forward 
time-contraction (~ O-SF^^Tj), and thus it masks the 
latter, as explained following eq. (IV. 10). This result of 



GW anti-beaming should be contrasted with the beam- 
ing of the EM signal, for which the forward signal is sig- 
nificant, and thus the radial time contraction sets the 
timescale. Fig. ([9| makes more apparent (for the same 
uniform jet openning angle and the same observer angle) 
the F-dependence of the effects: a larger F produces a 
stronger signal, which is also more strongly beamed for 
the forward observer (and thus reaches him sooner), and 
more time-contracted. 

A different acceleration mechanism - in particular, accel- 
eration driven by Foynting flux rather than a thermal 
fireball - would change the relation between F, r and t, 
and thus change the angular beaming and time contrac- 
tion characters of the signal. These should all produce 
different waveforms, with implications on detect ability. 
Thus from the GW signal we can learn about the accel- 
eration of the jet and about their angular distribution. 



9 



This could distinguish between uniform, structured and 
shotgun [45 jets (ah have different angular distributions) 
and possibly (after further study) between Poynting flux 
acceleration and thermal fireball acceleration. 
We notice an interesting feature of the wavforms observed 
from 6y < Oq (for example 0.05 in fig. shortly af- 

ter the signal h begins to rise, it drops sharply, reaching 
negative values, and only then it rises again. We observe 
this feature in a range of waveforms for various values of 
Oy^Oo (fig. (|8|). It is produced by the region of the jet 
surrounding the direction to the observer. In this region, 
the angle ^ around the FLOS produces phase factors e^*"^ 
(eq. (II.6)) which include negative contributions to 
Since this negative-dip arrives from directions close to 
the FLOS, it precedes the main signal. For sufficiently 
small observing angles, we find time periods where even 
the total signal /i+ becomes negative (fig. (|8|). A similar 
effect of phase-variation from angles close to the FLOS 
causes, for Oy ~ 6>o, the appearance of wiggles around the 
peak, as verified consistently for many waveforms and nu- 
merical simulations (fig. (|8|) 

By Fourier- transforming the waveforms from uniform jets 
(fig. (pT])), we find that the most prominent feature is a 
sine function, which matches our expectation of trans- 
forming a pulse in the time domain: the fourier com- 
ponents are approximately constant starting at very low 
frequencies and then they drop to zero. The width of 
the sinc^ which we measure as the first zero of the func- 
tion, corresponds to the width of the pulse in the time 
domain, and we find it at a frequency range /o ~ 600 Hz 
(for Tf = 30m^ and depending on the viewing an- 
gle, jet width and Lorentz factor). For a uniform jet 
of Oo = 0y = 0.1 we find (fig. ^) 




Mb) 






GWh(f} 




0=OA 




Observer eob,= 0.2 







FIG. 11. Fourier analysis of the waveforms produced by (a) 
uniform jets (F = 100, = 0.1) viewed from different angles 
Ov and (b) uniform jets (^o = 0.1, viewed from Ov = 0.2) 
for different F factors. Time is plotted in dimensionless units 

t/T/, and the fourier spectrum is scaled to h/{^^). All 
patterns match sine functions. 
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Uniform Jets, 60=0.1 



Y- 0.6783x0-7287 
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1000 



FIG. 12. F, the width (first zero) of the Fourier transform of 
/i, as a function of F, for Oy — Oq = 0.1. 



CONCLUSIONS 



/o 



0.67^ 



.0.73 



(IV.ll) 



We find that at larger viewing angles /o is independant of 
F and Oq , while for lower values of Oy it does not depent on 
Oy (fig. ([l3|). For Oy < Oq, we see a second (higher) range 
of frequencies for which the spectrum is non-negligible, 
matching the negative-dip in the waveforms for these an- 
gles. 



Fig.(14i) shows the waveforms for a structured jet (F = 
100,6^0 = 0.2) from several observation angles, display- 
ing similar features to a narrow uniform jet. It's fourier 
spectrum (fig.(14 3)) matches a sinc-pulse, regardless of 



the observation angle, indicating most of the signal arises 
from the very narrow uniform core. We also see that for 

Uy ~ "core ^ 

F ^, the signal includes a negative dip 
which follows the main signal rise, rather than precede it 
as in the uniform signature. 



^ As mentioned, in our model Tf in fact factors out linearly, and 
so all frequency results scale linearly with 1/Tf 



We studied the gravitational waves radiated from the 
accelerating ORB jets, examining both uniform jets and 
structured jets. We found that the quadrupole nature of 
GW causes anti-beaming away from the jet's axis, and 
thus the radiation peaks approximately at an angle F"-*^/^ 
off- axis. For a uniform jet, the signal is maximal just- 
outside the jet; for a structured jet, outside the core, but 
within the jet. For ^0 > F~^, the structured jet's GW 
signal is always more focused than its uniform counter- 
part, and as Oq increases, it remains more focused than 
the uniform jet signal (fig. ([3|, Q). As 7-rays from the 
jet are emitted within the jet width ^o, the GW peak is 
generally not within the GRB detectability cone; however 
the GW amplitude decreases only slightly away from its 
peak, and so a considerably high amplitude can be seen 
within the EM cone (fig. 

The waveform and the spectrum of the GW signal de- 
pend on the acceleration of the jet. We introduced an 
acceleration model of uniform linear acceleration, that 
is ba sed o n the Fireball model for thermal acceleration 
(eq. ( IV.4| )). It might be interesting to explore in the fu- 
ture other acceleration mechanisms and in particular the 
acceleration of a Poynting flux dominated jet. As the 
relations between F, r and t determine the beaming and 
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Sine width as a function of 6 




Uniform Jets, r=100 


. S(,= 0001 


90= 0.01 • 




X eo= 0.1 




♦ 90= 0.15 



Sine width as a function of Gobs ' 
Uniform Jets, 90=0.001 



0.0001 0.001 0.01 



0.0001 0.001 0.01 



TABLE L SNR for GW from GRB (FM = lO^^erg, at r 
bOOMpc) relative to expected noise in aLIGO & DECIGO 



FIG. 13. F, the width (first zero) of the Fourier transform of 
/i as a function of Ov, for (a) F = 100 and different jet widths; 
and (b) for different F's with width6'o = 0.001. 




FIG. 14. (a) /i as a function of time for a structured jet 
(F = 100, Oo = 0.2) at different observeration angles Oy. We 
see a peak pulse similar to a narrow uniform jet, and at very 
low angles Ov ^ Ocore ^ 

F we find the pulse is followed by 
a negative dip. (b) Fourier analysis of the same waveforms 
produced by structured jets (F = 100, = 0.2) viewed from 
different angles. The patterns match sine functions. 



time contractions, the different GW signals could teach 
us about the acceleration of the jet and could help dis- 
tinguish between Poynting flux acceleration and thermal 
fireball acceleration. 

For a source at distance r, of total energy M and a bulk 
Lorentz factor F, that accelerates over time Tj, the sig- 
nal's amplitude scales linearly with and the wave- 
form stretches over time linearly with Tf. The signal 
increases, first sharply over a timescale Uise ^ 0.5T~^Tf. 
this is fo llowed by a long slow increase until ~ 2Tf (figs 
( |6|7|8|9D ). The fast rise time represents a time contrac- 
tion of stemming from the angular quadrupolar anti- 
beaming of the GW, which stands in contrast with the 
F~^ time contraction factor of the EM radiation from 
the GRB. Another effect, we found, is the appearance 
of wiggles in the signal, preceding its peak (figs ( |8|9| )) , 
at observation angles within the jet cone, due to polar- 
ization. This implies that for a uniform jet the signal is 
much clearer and monotonous outside the jet's cone than 
the signal of a structured jet; however it also means that 
for any uniform jet detectable as a GRB (i.e., the ob- 
server is within the jet) we expect to see such wiggles in 
the GW signal, and they are more pronounced the closer 
the observe is to the jet axis. The waveform of a struc- 
tured jet displays a similar dip after the signal's peak. 
Thus the waveforms, and particularly the wiggles, offer 
insight into the internal structure of the jet, in particular 
to the angular structure and the acceleration. 
We turn now to analyze the detectability and exptected 
detection event rate using the frequency bands and sen- 
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sitivities of the aLIGO and DECIGO detectors. Ad- 
vanced LIGO 03 is planned to detect strain sensitivities 
of - 10-24 - 10-23/^;^^ in the range 30 - 300Hz, while 
DECIGO [471 SH] should reach the same sensitivity in 
the 0.1 — lOHz range. The characteristic frequencies of 
the GRB jet's GW's scale with 1/Tf (eq. (pVTl])). We 
have calculated the fourier signatures of different wave- 
forms, and found their main feature is a sine funtion 
(figs [n] [l2j [13) [m]). For a typical uniform jet of width 
60 = 0.1, F = 100 and acceleration time Tf = 30ms 
viewed at Oy = 0.1, we find the width of the sine is 
/o - 20HzT^-^^ - 600i^^. Following Flanagan et el [49]- 
53 , we have computed the Signal-to-Noise Ratio (SNR) 
for our waveforms using the expected sensitivity curves of 
aLIGO (with two operation-mode baselines: ZERO DET 
high P & ZERO DET low PgB]) and DECIGO (with two 
baselines: single cluster and 2-clutsers' correlation ^48]). 
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We used a typical parameter of total energy 
M = 10^^ er^, and considered a fiducial source at 
r = 500Mj9c, for comparison with the gravitaional 
radiation estimaes of [54]. Using the local GRB rate 
estimate of ~ 1.3Gpc~^i/r~^ [55]. our fiducial distance 
corresponds approximately to a yearly event. 
The results are shown in table [ij Using the baseline 
of DECIGO's 2-clusters' correlation and our typical 
parameters, we find SNRs in the range 5-100 (table [l|). 
As the SNR scales linearly with M and inversely with 
r, this indicates even sources 5-100 times more distant 
can be detected in this configuration of DECIGO. In 
particular, we should detect GW counterparts to most 
of detectable GRBs. We also expect to see such GW 
signals without seeing the corresponding GRB itself, 
due to their different beaming properties: the GW's 
amplitude is considerable outside a ~ r~^/^ cone, ~ 50 
times wider than the 7-ray angl^ For the signals 
coming from jets pointing towards us at angles that are 
of order r~^/^ we expect significant orphan afterglow 
[54] signals to accompany the GW signals. 
A GRB jet's GW is also expected to be preceded by 
GW from its progenitor, and the coincidence of the two 
signals (even on different frequency scales and using 
different detectors) could improve the chances of finding 
them. 

Using aLIGO or a single DECIGO cluster and our 
typical parameters, the SNR for a uniform jet is ex- 
pected to be ~ 1 — 5 • 10~^, while for a structured 



jet it can reach 10~^; both are undetectable. A GW 
signal can be detected together with a GRB if it is ~ 10 
times stronger /closer, thus ~ 1000 times less frequent 
(~ 10~^yr~^). This rate is prohibitively low even 
considering isotropic GW signals (whose GRB and/or 
orphan afterglow counterparts are not visible). 
Other astrophysical phenomena might produce similar 
GW signals. In particular. Micro Quasars ejecting mass 
and energy at relativistic velocities are less powerful 
{TM - lO'^^erg), but they are much closer to us 
(r ~ 5Kpc). Using eq. (II.7), we estimate their peak to 
be h ~ 10~^^, two orders of magnitude less than the 
typical GRB signals; we expect their typical frequencies 
to be / < 35Hz (for a similar Tf but much smaller 
r ~ 2, compare eq. (|IV.llD,fig. (pt). Active Galactic 
Nuclei can also display acceleration (magnetic) and 
emit much more energy, over longer time scales and 
lower frequencies. The same techniques, once adapted 
to continous rather than impulsive sources, could be 
applied to analyze the GW output of such sources. 
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Appendix A: Additional Resources 

All of the code used is available at 
http: / /www. phys.huji.ac.il/~ofek/GWGRB/ 
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